Fluorescence measurements of expanding strongly-coupled neutral plasmas 
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We report new detailed density profile measurements in expanding strongly-coupled neutral plas- 
mas. Using laser-induced fluorescence techniques, we determine plasma densities in the range of 10 
to 10 9 cm -3 with a time resolution limit as small as 7 ns. Strong-coupling in the plasma ions is in- 
ferred directly from the fluorescence signals. Evidence for strong-coupling at late times is presented, 
confirming a recent theoretical result. 
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Strongly-coupled Coulomb systems appear in extreme 
conditions, such as in quark-gluon plasmas 0, in laser- 
fusion plasmas 0, and in some astrophysical settings. A 
new class of strongly-interacting neutral plasmas was re- 
cently demonstrated using the tools of laser-cooling and 
trapping 

SUSIE II These "ultracold" neutral plasmas 
occupy a unique position in phase space. In these plas- 
mas it is possible to create strongly-interacting Coulomb 
systems at modest densities because the initial electron 
and ion temperatures can be in the milliKelvin range. 
The initial ion-ion and ion-electron interaction strength 
can also be selected with great precision. 

Recent experimental work in this field has used absorp- 
tion imaging techniques to make temperature and den- 
sity measurements in expanding ultracold neutral plas- 
This work explored the 50 to 1000 ns time pe- 



riod after plasma formation in great detail. Correlation- 
induced heating was observed in the plasma ions. The 
ion coupling parameter, given as the ratio of nearest- 
neighbor Coulomb energy to the kinetic energy, equili- 
brated just inside the strongly-coupled regime, with a 
coupling parameter around 2. 

Radio-frequency (RF) excitation techniques have also 
been used to determine the average ion density and the 
electron temperature in these systems 0, H, |tj. These 
studies confirm theoretical predictions regarding the gen- 
erally self-similar Gaussian expansion of the ions and 
the clamping of the electron temperature in the weakly- 
coupled regime. 

In this letter we report laser-induced-fluorescence mea- 
surements of ions in expanding strongly-coupled plasmas 
as a tool to study the spatial and temporal evolution of 
the ion temperature and density. This measurement tech- 
nique has a 7 ns temporal resolution limit. We measure 
plasma densities as low as 10 5 cm" 3 at effective plasma 
temperatures of 100 K. The maximum density that can 
be measured is limited by radiation trapping, and for 
spherically-symmetric systems in the milliKelvin range 
is limited to around 10 9 cm" 3 . The temporal resolution 



and dynamic range of this method in ultracold plasma 
measurements surpass those currently seen in absorp- 
tion spectroscopy, and rival the sensitivity of RF spec- 
troscopic methods. 

Much of the experimental setup has been described 
previously ^(j- We create a calcium magneto-optical 
trap (MOT) using the resonance transition at 423 nm. 
Up to 50 mW of 423 nm radiation is generated by 
frequency-doubling a diode laser system at 846 nm using 
a periodically-poled KTP crysal in a resonant build-up 
cavity. The 423 nm MOT light is detuned one line-width 
below the atomic transition, making the MOT tempera- 
ture around 1 mK. The density distribution is approxi- 
mately Gaussian in x, y, and z, of the form 



= n cxp[~(a: 2 + y 2 )/al - z 2 / 1 ] 



(1) 
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The 1/e 2 radius is oq = 0.5 mm in the x and y dimen- 
sions, and Pq = 2.5<7o = 1.25 mm in the z dimension. 
We use a grating-stabilized 657 nm diode laser as a re- 
pumper to increase the number of atoms in the trap. The 
peak density determined by absorption measurements on 
the atoms is 4 x 10 9 cm" 3 . 

We photo-ionize the atoms in the MOT using a two- 
color, two-photon ionization process. A portion of the 
846 nm diode laser radiation from the MOT laser is 
pulse-amplified in a pair of YAG-pumped dye cells and 
frequency doubled. This produces a 3 ns-duration laser 
pulse at 423 nm with a pulse energy around 1 /zJ. This 
laser pulse passes through the MOT, and its peak in- 
tensity is a few thousand times greater than the satura- 
tion intensity. A second YAG-pumped dye laser at 390 
nm counter-propagates the 423 nm pulse and excites the 
MOT atoms to low-energy states above the ionization 
potential. We photo-ionize 85-90% of the ground-state 
atoms in the MOT. The minimum electron temperature 
is limited by the 30 GHz (measured) bandwidth of the 
390 nm laser to about 1 K. 

Ions in the plasma scatter light from a probe laser 
beam tuned to the Ca II 2 Si/2 — 2 P\ii transition at 397 
nm. The probe laser is generated by a grating-stabilized 
violet diode laser locked to the calcium ion transition 
using the DAVLL technique m a large-bore, low- 
pressure hollow cathode discharge of our own design. The 
probe laser is spatially filtered and a few /zW of power 
is focused to a Gaussian waist of 130 /im in the MOT. 
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FIG. 1: Laser-induced-fluorescence of ions in an expanding 
ultracold plasma. In the upper panel, the probe beam is par- 
allel to (and collinear with) the long axis of the initial plasma 
density distribution. In each successive panel, the probe beam 
is moved 0.22 mm = 0.44<ro farther away from this symmetry 
axis, where <ro is the 1/e 2 transverse Gaussian radius of the 
initial plasma density distribution in the x dimension. Note 
that the time scale is logarithmic and the fluorescence scale 
is linear. The solid line is the expected signal from a thermal 
expansion model, as explained in the text. 



a steady excited-state population after a few radiative 
lifetimes. We also see this same turn-on time when we 
perform absorption measurements on the plasma ions. 
A higher probe laser intensity would drive the ion pop- 
ulation into the excited state more quickly, making it 
possible to probe the plasma at earlier times. 

After the signal peaks, it rapidly decays. At t ~ 100 ns 
the decay slows down. This is due to correlation-induced 
heating in the ions, similar to previously published ab- 
sorption measurements 7] . During the first microsecond 
of the decay, while the spatial distribution of the plasma 
ions has not changed, it is straightforward to convert 
this fluorescence decay signal into velocity as a function 
of time. With a few assumptions, we can determine the 
temperature, density, and strong-coupling parameter of 
the plasma ions. 

An atom moving with velocity v z parallel to the 
propagation direction of the probe laser beam has a 
Lorentzian probability of scattering a photon, propor- 
tional to l(v z /b) 2 + The constant b — 7A/2 is the 
velocity that corresponds to a Doppler shift equal to the 
natural linewidth, 7 = l/27rr is the natural linewidth of 
the transition, and A is the transition wavelength. The 
initial velocity distribution is Maxwellian, of the form 
exp(— v1/2vj h ) / y/irvth- For systems in thermal equilib- 
rium, vth is the thermal velocity. As mentioned in the 
literature 0, the temperature should be a measure of 
the random motion of ions. In these kinds of measure- 
ments for expanding ultracold plasmas, this random mo- 
tion is observable only at early times. The random mo- 
tion is quickly overwhelmed by the directed motion of 
the plasma expansion. So Vth quickly looses its meaning 
in a thermodynamic sense. However, when the ions are 
"heated" by correlations the Maxwellian approximation 
is still valid. 

We integrate the product of the Lorentzian line shape 
and the Maxwellian distribution to get an expression for 
how the fluorescence signal, s(t), changes with the width 
of the ion velocity distribution as a function of time: 



We average repeated measurements of the scattered 397 
nm radiation with the probe laser in a given position. 
This produces a time-resolved signal proportional to the 
number of atoms resonant with the probe beam in a par- 
ticular column of the plasma. By translating a mirror 
just outside the MOT chamber, we scan the probe laser 
across the ion cloud. In this manner we obtain tempo- 
ral and spatial information about the plasma expansion. 
Typical fluorescence measurements are shown in Fig. ^ 

When the probe beam is on axis (Fig. ^ top panel), 
the fluorescence signal takes roughly 20 ns to reach its 
maximum value. This duration is a few times the 7 ns ra- 
diative lifetime of the Ca II 2 P\/2 level. This is an exam- 
ple of the classic Rabi two- level atom problem with spon- 
taneous emission, where a collection of atoms driven on 
resonance well below the saturation intensity will reach 



s(t) oc — exp fi-] erfc ( — | (2) 
«th WJ \V2vth) 

where erfc(x) is the complimentary error function of the 
parameter x. If v t h is known at a particular time, this 
equation can be inverted to give the change in the ther- 
mal velocity as a function of time. The initial neutral 
atom velocity is near the Doppler limit at 40 cm/s. The 
neutral atoms also experience an ionization recoil veloc- 
ity of 40 cm/s. In the inversion of Eq. |2 we set the 
scale so that a linear extrapolation of the signal to time 
t = gives the quadrature sum of the Doppler and recoil 
velocities, 56 cm/s, as shown in Fig. |21 

A thermal model of ultracold plasma expansions was 
presented in Ref . [l2l fl3j . We have analytically solved 
the coupled differential equations presented in that work 
and find the analytical solution to the directed ex- 
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FIG. 2: The mean ^-component of the ion velocity as a func- 
tion of time after plasma formation. The points show the 
velocity extracted from the ion fluorescence signal shown in 
the top panel of Fig. The solid line is a fit of the data to a 
thermal model [T^l . as described in the text. The dashed line 
shows the extrapolation to 56 cm/s at time t = 0. The inset 
shows the ratio of the data to the thermal model, accentuating 
the correlation-induced heating signal. 

pansion velocity to be 



where a(t) = \/ <J 2 + v 2 t 2 is the time-dependent width of 
a self-similar Gaussian expansion. At early times, when 
v 2 t 2 <C cr 2 , the velocity reduces to v(r, t) w v 2 t r/a 2 . 

In applying this expression for the expansion velocity 
to the data in Fig. |2 some care is required. The model 
assumes ions initially at rest. However, in our experi- 
ment the plasma ions have an initial electrical potential 
energy which is quickly converted to kinetic energy (see 
the discussion of correlation heating below). This gives 
an offset to the t = velocity predicted by the model. 

Otherwise, the application is straightforward. The 
probe laser interacts with ions in a particular volume. 
During the first microsecond, the large-scale ion distri- 
bution does not change, even though the ions are ac- 
celerating. The expression in Eq. |3| predicts a velocity 
that depends linearly in time, with v 2 r/a 2 averaged over 
the part of the density distribution interacting with the 
probe laser beam. A fit of our data to this model is 
shown in Fig [21 For the early expansion times, the lack 
of spherical symmetry in our plasma is unimportant in 
fitting this data, because the average v 2 r/a 2 appears as 
a fit parameter. 

This figure illustrates correlation-induced heating, and 
provides an independent check of the plasma density. 
When the plasma is created, the ions quickly move from 
their initial disordered state into a state with greater or- 
der. This motion is initially synchronized at the plasma 



frequency. But because of variations in the local elec- 
tric field at each ion, the oscillatory motion becomes un- 
synchronized. The time from plasma creation to the max- 
imum temperature (minimum signal) in this oscillation 
is equal to t — l/4w p , where uj p is the average plasma fre- 
quency. As seen in the inset to Fig. the peak occurs at 
60 ns, giving an average plasma density of 4.0 x 10 s cm~ 3 , 
about a factor of three lower than what we expect based 
on the ionization fraction of the neutral MOT. The cause 
of this discrepancy is not readily apparent. Our numeri- 
cal modeling of the plasma oscillation dynamics indicates 
that the average density determined from the observed 
oscillation period reproduces the average plasma density 
to within 20%. Perhaps this discrepancy indicates the 
number of plasma ions that have recombined into high- 
lying Rydberg states at early times. 

The strong-coupling constant in this system can be de- 
termined from data extracted from Fig. [S] The kinetic 
energy of the system is equal to mv 2 /2, where v 2 = v 2 + 
Vy + v 2 , and the random thermal motion due to correla- 
tion heating is assumed to be isotropic. The temperature 
is determined from 3fcfcT/2 = mv 2 /2. Using the derived 
values of v z — 8 m/s and n = 4 x 10 8 cm -3 , the strong- 
coupling parameter is T = (e 2 /47reofcsT)(47rn/3) 1 / 3 = 4. 

During the time from 1 to 40 /zs, the fluorescence signal 
depends on the time-evolution of the ion temperature, 
the overall plasma expansion velocity, and the density. 
As the plasma expands, the number of ions in the probe 
laser beam decreases. Also, as the velocity distribution 
widens, the fraction of atoms in the beam that are in 
resonance with the probe laser frequency also decreases. 

The thermal expansion model in Ref. |13| provides a 
theoretical framework for interpreting our data. Most 
of the ions in the plasma lie in a region in which the 
electric field depends linearly on the spatial coordinate. 
The ansatz presented in Ref. ^2] simplifies the system 
by imposing this linear dependence of the electric field for 
the entire system. While this is clearly not true for ions 
near the edge of the plasma, the number of ions in that 
region is small. Sophisticated approaches to modeling 
this system have shown that the thermal model captures 
the major features of the expansion dynamics 

The lack of spherical symmetry in our experiments 
changes the slope of the electric field in the z direc- 
tion relative to x (or y), because the plasma distribu- 
tion is elongated in z by a factor of 2.5 (see Eq. QJ. 
Keeping with the ansatz of Ref. 0], we assume that 
the Gaussian plasma expands in a self-similar fashion. 
Similar to Eq. the expansion velocity in the x di- 
rection is v x (t) = v 2 t xja 2 {€). The expansion veloc- 
ity in the z direction is v z (t) = v 2 t z/P 2 (t), where 
(3 = y/ (2.5(7o) 2 + v 2 t 2 . This result, which relates the 
velocity to the position, will be important in modeling 
the fluorescence signals, as discussed below. 

In the experiment, the probe laser beam is spatially fil- 
tered and focused into the plasma with a Gaussian beam 
waist of w — 130 fim. The position offset of this probe 
laser relative to the plasma is denoted by the distance a. 
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The two-dimensional convolution of the probe laser and 
the plasma distribution is proportional to the number of 
plasma ions per unit length along the probe beam, 

v = p( W "+ 2g») exp (~^w) cxp irw) • (4) 

Multiplying this expression by the Lorentzian lineshape 
gives a probability that an ion in the probe beam column 
moving with a z-componcnt of velocity v z interacts with 
the probe laser. We substitute for v z the expression de- 
rived previously, and integrate Eq. 0]over z to derive the 
predicted time-dependence of the fluorescence signal: 

, (l\ wH ( 2a 2 £ 2 \ 

* (t) " 6rfC U J 0(w 2 +2a 2 ) ° XP (f 5*+2? + J 2 ) 

(5) 

where I = 'yXa 2 (t)/2tv 2 . This derivation neglects the in- 
fluence of the thermal velocity and is valid at times where 
the directed expansion velocity dominates the thermal 
velocity (a regime opposite to that for Eq. 0). Adding 
the thermal velocity to the model would further reduce 
the predicted signal levels at early times. 

This fluorescence model is the solid line plotted with 
the data in Fig. ^ The only two adjustable parameters 
are the expansion velocity v e and the overall amplitude. 
They are chosen so that the model agrees with the data 
in the last panel of the figure. Considering the simplicity 
of this model, it is surprising that it matches the general 
shape and order of magnitude of the fluorescence signal. 

Relative to the model, the experimental data are 
consistently higher at later times. Additionally, data 
recorded with the probe laser beam between a = <jq 
and 3(7o show structure that cannot be explained by the 
model. Predictions of ion-acoustic waves freezing into 
the plasma expansion have been made. However, the 
observed modulations are not stationary relative to the 
plasma expansion. A vertical slice through all of the 



panels in Fig. Ogives the density profile at a given time. 
The data show that two local maxima that develop in the 
interior of the plasma coalesce into one another and dis- 
appear relative to the smooth Gaussian shape predicted 
by the model (lowest frame). Contrary to predictions, we 
see no evidence for shock formation in these plasmas. 

It is likely that the plasma ions are strongly coupled 
throughout the expansion. The ions equilibrate after 100 
ns with a coupling parameter of T = 4. Recent hy- 
brid molecular dynamics calculations |15| show that af- 
ter an initial oscillation, the coupling parameter always 
increases. The absence of a shock front in the plasma 
expansion and also the observed suppression of local den- 
sity variations is evidence that the plasma ion system is 
collisionally stiff, once again suggesting strong coupling. 

In conclusion, we have demonstrated a new technique 
for determining the the density of ions in an expanding 
ultracold plasma using fluorescence spectroscopy. Our 
implementation of this technique has a time resolution 
of 7 ns. We report measurements of the density pro- 
file for up to 50 /is, and demonstrate a sensitivity limit 
around 10 5 cm~ 3 . This method surpasses the time res- 
olution and dynamic range of previously reported tech- 
niques. It rivals the sensitivity limit of RF techniques, 
and provides detailed spatial information of the plasma 
density profile. We have created a plasma just inside 
the strongly-coupled regime, and observed changes in the 
plasma density over time. Based on recent theoretical 
work, it is likely that these plasmas are strongly-coupled 
for the entire measured expansion. The general features 
of our measurements are in reasonable agreement with 
a thermal model, although the remaining discrepancies 
require more sophisticated methods to interpret quanti- 
tatively. 
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